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We solve the problem of A'' non magnetic impurities in the staggered flux phase of the Heisenberg 
model which we assume to be a good mean-field approximation for the spin-gap phase of the cuprates. 
The density of states is evaluated exactly in the unitary scattering limit and is proportional to 
1/ (ujln^{\uj\/ D)"^ , in analogy with the ID case of doped spin-Peierls and two-leg ladder compounds. 
We argue that the system exhibits a quasi long-range order at T = with instantaneous spin-spin 
correlations decreasing as rii/ln^ (riiRij) for large distances Rij and we predict enhanced low energy 
fluctuations in Neutron Scattering. 



It is now established that the normal phase of the un- 
derdoped cuprates possesses a magnetic gap with the 
same anisotropy as the d-wave superconducting gap. The 
essential physics of the high-Tc superconductors can be 
captured by focusing on the Cu02 planes. A microscopic 
starting point of theoretical analysis is the so-called t-J 
model. In a slave-boson mean-field formulation, the sys- 
tem undergoes spin and charge separation: an electron 
in these highly correlated materials is a composite object 
made of a spin 1/2 neutral fermion (spinon) and a spinless 
charged boson (holon). The gap in the magnetic excita- 
tions of the normal state can be viewed as singlet for- 
mation between pairs of neutral fermions in the absence 
of coherence between the holons. We assume that this 
mean- field picture captures essentially the initial idea of 
Anderson 0| of a Resonance Valence Bond ground state 
for the normal phase of the cuprates. 

Substitution of Cu ions in the conduction planes of 
high-Tc cuprates by different non magnetic ions presents 
an important experimental tool for the study of the 
metallic state. Unusual effects has been revealed espe- 
cially when the materials were doped with Zn. The va- 
lency of Zn is Zn^+ and compared with the Cu^+ 
case, one electron is trapped by an additional positive 
charge of the nucleus, forming a singlet at the Zn site. 
In the spin-gap phase, it is experimentally found by 
NMR experiments that a local magnetic moment 
of spin 1/2 appears on Cu sites neighbouring the Zn im- 
purity. 

In high-Tc cuprates the persistance of AF fluctuations 
in the metallic state is probably one of the most striking 
feature. Inelastic Neutron Scattering have established 
previously the existence of an energy gap in the imagi- 
nary part of the dynamical susceptibility in the normal 
and superconducting phases of the pure compound |^] 
(without Zn). Substitution of Zn in Cu02 planes shows 
a striking transfer of spectral weight from high to low 
energies, partially filling the spin-gap This is the 

signature of strong enhancement of AF correlations in 
the spin-gap phase. 

On the theoretical side, impurity-induced moments 



have been studied for a variety of quantum disordered 
spin systems. Especially in one dimension, doped 
spin ladders P, p^ , pT 12 1 such as SrCu2 03 and spin- 
Peierls 0,11^ systems such as CuGeOa have been in- 
vestigated. These systems have in common with the un- 
derdoped cuprates the presence of a gap in the spin exci- 
tation spectrum. It has been shown theoretically as well 
as experimentally jl^ that Zn doping is also the origin of 
a strong enhancement of AF correlations in the ground 
state. The picture is that the impurity breaks up the sin- 
glet, leaving behind a 5 = 1/2 local moment. Moments 
on the same and opposite sublattices interact with fer- 
romagnetic and antiferromagnetic exchange respectively. 
There is no frustration in this system. The problem has 
been mapped onto the random S' = 1/2 chain with ran- 
dom sign interaction Q and is believed to exhibit long 
range or quasi-long range order at T = 0. It has been 
suggested that similar physics may operate in two dimen- 
sional doped t-J model |]ll|] . 

Earlier work has treated the effect of a single non 
magnetic impurity in the 7r-flux state |]l5| and the pair- 
wise interaction between them [ p^ . It was found that 
each impurity creates a bound state in the pseudogap 
at = 0. Interaction between a pair of impurities 
leads to a level splitting between these states given by 
A(i?) - l/(i?ln(i?)). 

Here we solve the problem of N non magnetic impuri- 
ties in the staggered flux phase of the Heisenberg model. 
We find that the overlap between the bound states leads 
to a broadening of the (5-function peaks and the density 
of states is given exactly by 

5p{u) = n,/ {\Lo\\n^\uj\lD)) , 

{rii being the density of impurities). We also confirm that 
because the staggered flux phase introduces no frustra- 
tion between the local moments, a quasi-long range stag- 
gered order exists at zero temperature. We present our 
results in terms of a 7r-flux phase even though our con- 
clusions are general for any staggered-flux phase. In our 
view the quasi long range order found for the 7r-flux phase 
at half filling accounts for the general enhancement of an- 
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tiferromagnetic correlations found by INS experiments in 
the underdoped compounds. 

In the absence of impurities, the Hubbard model 



H — —t ^ (c|cj + h.c.) + U ^ ni-^riii 

(id), a- i 



(1) 



for I] jt 1 at half-filling may be canonically trans- 
formed and projected into the nondoubly occupied sub- 
space to become the S —\ AF Heisenberg model 



AF 



(2) 



where J = At^ /U . The mean-field approach assumes a 
nonzero expectation value for Xij — (clo-Cjo-}, which is a 
link variable of a U(l) lattice gauge field coupled to the 
remaining fermions Xij — exp[i A.dl]. We parametrize 
the TT-flux phase by taking Xij = xe*^'^ i where Oij = —Oji 
and 9ij = if ij || x, and ±-| if ij || y. 

Dividing the square lattice into two sublattices 
A and B, the mean- field hamiltonian can be di- 
agonalized via the canonical transformation ~ 

'''{ipk,o - ipk, i), with 
3^ + cos^ ky , leading 



ipk,o + "0^,1 and = e 

giVfc — (cos fca; -I- i cos fcj,) /-^/c 

to Hmf = T.keBZ,a=o,i^k{-'^)°'^l^a^k,a, whcrc the 
summation over k extends over the full Brillouin zone 
and the constant terms have been omitted; ±i5fe — 
± ^J cos^ kx + cos^ ky are the quasiparticle eigenvalues 
that will be linearized around the four nodes (± J, ±^). 

The impurities are treated as repulsive scalar poten- 
tials randomly distributed on the lattice 



N 



imp 



(3) 



We assume that the unitary limit {Vq 1) is physical 
because the electrons on Cu sites are really trapped by 
Zn impurities in experimental realizations. 

We define the Green's function G^^, (iuJn) = 
-{T'rV'fc,a(T)V'I,_a,(0))c^„, whcrc the brakets ( de- 
note the Fourier transform in Matsubara frequencies. 
Note that no impurity averaging has been made. The 
bare propagators for quasiparticles satisfy ^{iujn) = 



lUJr,-{-l}''Ek • 

The equation of motion for the Green's function can be 
closed algebraically: it involves N unknown propagators 
(corresponding to the scattering by each impurities) that 
can be evaluated self consistently. This leads us to the 
following T-matrix equation 
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G\.a^k,k'5a,, 



(4) 



QlJ-B.,)T,j{iio^)gl.A-B.,) 



where G^^ai^i) = Gfc,Q ^xp (ik • R^) if Rj belongs to 
sublattice A and 5" ^(R^) = ^ exp (ik • R^ — iipk) if 
Ri belongs to sublattice B. The T-matrix is given by 
T = —VqM^^, with M the matrix that closes the equa- 
tions of motion. M is a. NxN matrix {N is the number of 
impurities) each entry of which corresponds to the scat- 
tering from one impurity to another. It can be written 
by blocks as 



M = 



An/2 



(5) 



where A and D describe the scattering by impurities in 
the same sublattice and B the scattering by impurities in 
different sublattices. We have noted in indices the rank of 
the matrices which is roughly N /2 but the solution works 
for any macroscopic number of impurities in each sublat- 
tice. Dividing again in the same way each sublattice into 
two sub-sublattices Ai, A2 and Bi, B2, we define the co- 
efficients Aij — Sij + Vo A{iuJn, R-jj); if R-i and Rj belong 
to the same sub-sublattice and Aij — when R; and Rj 
belong to different sub-sublattices. Bij — Vq B{iLUn,'R,ij) 
and Dij = Aij. Two impurities in the same sub- 
lattice are related by the propagator A{iujn , Rjj ) = 
Jf J2k,n Gfc,n(*'^n) ^xp (ik • Kij) , with Af the number of 
sites in the lattice. When two impurities are in dif- 
ferent sublattices, they are related by S(jw„,Ry) = 
7^Efe,nG°,„(jw„)(-l)"exp(-i(^fc-|-ik-R,j). Due to 
the symmetries of the 7r-flux phase, there is no interac- 
tion between two sub-sublattices that belong to the same 
sublattice. 

The coefficients A{iLUn, Ry ) and B{iujn, Rjj) are given 

by 



yl^(iw„, R) 



R\^n\ 

D 



where $"^(i?) = ^ (cos(7r_Rj;) -I- 1) and Kq is the modified 
Bessel function of rank zero; 

6(iw„, R) ~ 4i exp [iQs • R + (i?) ^i^i ^''^^ 



D 



with 

Q2 = (f,f), $^(i?) = i(l-cos{7ri?,))exp(-2i(^), >f 
being the angle between R and x . 

In the regime where Rij\uJn\ ^ D, Aij{iLOm'R) ~ 
4/(7r_D'^) In |i?a;„/Z?|. This logarithmic structure is 
caracteristic of the system of Dirac fermions in 2 D. It 
is the origin of the breakdown of the usual perturbative 
expansions jT^ ]: logarithmic divergences will indeed ap- 
pear in all orders in the diagrammatic expansion, making 
it necessary to re-sum the whole series of diagrams in or- 
der to control any calculation. 

Here we take a non perturbative approach and directly 
invert the matrix M. We notice that Ajv/2 approaches 
unity in the limit of low frequencies. We use 
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^B{iuJn, Rik)B{iuJn, Rkj) 



ln(|c^„|i?,,/I?) S , (6) 



where S" is a N/2 x N/2 matrix satisfying Sij = 1 if Rij < 
D/\ujn\ and elsewhere. The matrix product *BB 1. 
After taking the limit Vb — > oo, we invert the matrix M 
by blocks: 



T{iLU„) = 











(7) 



Impurities interact only if they are on different sublat- 
tices. We have supposed that B is a square matrix 
but the proof can easily be generalized for a non square 
one 

We proceed to evaluate the additional density of 
states Sp{uj) = — ^ -f'^i^j, ^ (JG^^ (w + z<5), where 6G 
is the part of the Green's function in eq. (0) due 
to impurities. It can be read from eq. (m that 

5p{Lj) = -i Im Tr{Td/dujM^, so that 5p{uo) = 

— |: Im d/duj{\nDetB). The idea is now to evaluate 
this trace without exactly calculating the T-matrix (H). 
We know that DetB = 1/2 Det{*BB) and we use 
the logarithmic form of *_Bi? to get DetB ~ 
1/2 f{\Lj\) {\uj\/D) where / is a prefactor which 
doesn't depend on the number TV of impurities. The 
logarithmic structure of the problem that invalidates 
the usual perturbative expansions has been used here 
to factorize the leading divergence in the evaluation of 
Det{*^ Bli). Since the prefactor f{\Ld\) doesn't depend on 
iV, its contribution to the density of states will be neg- 
ligible in the thermodynamic limit. Hence in the ther- 
modynamic limit and after analytic continuation, we get 
5p{ijj) — -^'^^ (in In |-^|) and per unit of volume, 
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2\uj\ \n^{\uj\/D) + [T:/2f 



(8) 



Note that 6p{uj)duj — ni/2, there is thus as many 
states created at zero energy as impurities in the system. 
For one impurity alone jlsllBl, a (5-like bound state is cre- 



ated at w = 0. If the impurities were totally uncorrelated, 
we would find N ^-functions at a; = 0. Here we see clearly 
the overlap of impurity states which leads to a broaden- 
ing of the 5 functions. It is interesting to compare this 
density of states with the one found for one dimensional 
spin-gap systems such as spin-Peierls Jl3| and two-leg lad- 
ders 0Jll|,§ systems where 5p{uj) - 1/ {\lo\ ln^(|u;|/D)). 
It is also worth noticing that the very same form of the 
density of states as eq. (^) has been obtained in other 
apparently unrelated problem p9| ] for reasons which are 
not clear to us so far. 

In our method, we in fact calculate exactly many ther- 
modynamic quantities. For instance the uniform mag- 
netic susceptibility is given by 



X(T) = /3 / de5p{e) 



1 



2cosh'(f ) 



(9) 



From the asymptotic behavior of 5p{e) at small e we 
find that at low temperatures x{T) ~ n,/ (2Tln(l/T)), 
which can be interpreted as a Curie-like behavior with a 
vanishing Curie constant C{T) = rii/ (21n(l/r)). 

The free-energy 



AF{T) 



-2T 



deSp{e) ln(l 



leads to a specific heat CviT) ^ l/(ln^(l/T)) and a low 
temperature entropy S{T) ~ l/ln(l/T). However the 
entropy vanishes very slowly with temperature. 

We address now the important question of long range 
order in the system by computing the instantaneous cor- 
relations {S+S~) = lim^^o+ {Sti'r)Sy{0)). It is easy 
to show that (S^S^) ~ lim^^o+ (G'Jj(T)G]j(-T)), where 
Gij (r) is the Fourier transform in real space and com- 
plex time of the exact Green's function of eq. (||). We 
verify that this two-spin correlation is staggered. If 
i and j belong to the same sublattice, dGij{iujn) — 
X]M^(jWn,Rjfc) Tki B{iLUn,'Ri]) which owing to sim- 
ply renormalizes the bare Green's function: 5Gij{iuJn) = 
—niA{iuJn, Rij)- When i and j belong to different sub- 
lattices, there is one non trivial quantity SGij{iuJn) = 
J2ki-^i'i-^n,'Rik) Tki A{iuJn,'R.ij)- A simplification is ob- 
tained after noticing that A = B^^ {dB /diojn), which 
leads to 



(10) 



(5Gy(iw„) = [dB/diLUn] [dB-ydiuj. 



Writing dB^/diuin = 2B {dB/diuJn) and evaluating 
this product with the use of eq. (|6|), we show that 
dli"^ /diun ^ {^\T^{\uJn\/ D)/diuJn) S. Here again we have 
factorized the logarithmic divergence in B^ and shown 
that the derivative with respect to iuJn acts mainly on 
this factorized logarithm. Consequently dB~^ / diojn ~ 
{d{l/ hi\LOn\/ D) / diuJn) S^^. We now make an approxi- 
mation for the inverse of S. We take S~j^ = {\uJn\/D) Uij. 
For all j inside a circle of radius D/\uJn\ around the 
point i, Uij is a random configuration of ±1 so that 
D/\ujn\- In addition all the points j situated 
in the external boundary of this circle have Uij = ~l/n. 
Elsewhere Uij — 0. The main difficulty in the inversion 
of S is that two circles of radius D/\uJn\ centered around 
two points i and j very close to each other will overlap, 
leading to the same number of non zero coefficients in 
the lines i and j of S. In order to differentiate the sums 
X)fc '^'ife'S'^r/ and '^k^ikS'/I^ we have thus used the ex- 
ternal boundary of the circle to compensate its volume 
in . Using this inverse we perform the evaluation of 
eq. (10). In order to take into account the vertex correc- 



tions due to scattering processes, we perform the average 
over disorder after having evaluated the product of the 
two Green's functions. We find that at T = 0, 
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ists-) 



(In^ + (^/2)2) 



(11) 



for impurities in different sublattices. N non magnetic 
impurities randomly distributed in a 7r-flux phase stabi- 
lize a quasi long range staggered order. This result has 
to be compared with the pure 7r-flux instantaneous cor- 
relations which decay algebraically in l/i?"* |Q at zero 
temperature. Physically we understand what happens by 
introducing impurities one after another. For one impu- 
rity, a local moment with S = 1/2 is created in the vicin- 
ity I p^ . Under on-site repulsion, two impurities interact 
via an effective Heisenberg exchange H = J(i?)Si.S2, 
with J{R) ^ l/{R\n{R)) if they are located on different 
sublattices. However two impurities located on the same 
sublattice don't interact In our model, eventhough 
the on-site repulsion is treated at the mean-field level, the 
one particle wave-function remains commensurate with 
sublattices. This supports the intuition that no frustra- 
tion will be introduced in the model when a macroscopic 
number of impurities interact with each other. This ef- 
fect is the origin of the quasi long range order that we 
find. 

A similar picture is obtained for the ladder or spin- 
Peierls cases except that in these cases the instanta- 
neous spin correlation decreases exponentially with den- 
sity, reflecting the fact that interactions between local 
moments fall exponentially with the distance. Here 
the instantaneous correlations on different sublattices 
are proportional to the density of impurities because 
a pair of impurities interact via an effective exchange 
J{R) ^ l/(i?ln(i?)). We suggest that the logarithmic 
decay in ( [Til ) is due to quantum fluctuations at zero tem- 
perature even though the precise power in the logarithm 
may depend on our approximation. 

What would be observed if a Neutron Scattering ex- 
periment would be performed on this system? We 
have computed here the structure factor S{Q, f2) = 
l/-^X]g x"(Qj ^) integrated around (7r,7r). For two 
points i and j in different sublattices, Xij(Qi^) 

pressed as Xij(Qj^^ + *^) — l/l^r/?) Jq duj ImGij{uj — + 
iS) ImGji [ijj -I- iS) . Using the same tricks as before for 
the instantaneous correlation function, we show that at 
low frequencies 



\n/D\ 



(12) 



Even though we don't have true long range order in the 
system (which would have led to S'(Q, 17) = 5{ri) ), Neu- 
tron Scattering experiments would see a divergence of 
the integrated staggered structure factor at low frequen- 
cies, signaling a strong enhancement of antiferromagnetic 
correlations. 

The situation will change for experimental systems 
with finite oxygen doping. If we take as a theoretical 



starting point the d-RVB ground state, the effect of oxy- 
gen doping will be to move the nodes away from the four 
points (±7r/2, ±7r/2) in the canonical spectrum. Excita- 
tions from node to node will then slightly differ from the 
commensurated vectors (tTjTt). We believe that instead 
of quasi-long range AF order, this mechanism may give 
rise to an incommensurate quasi- long range order. Fur- 
thermore, it is likely that the finite frequency response 
integrated over g-space will not be that different from 



that given by eq.(12), which may explain the enhance- 



ment of AF correlations observed in Zn doped YBaCuO 
by neutron scattering 
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